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Statement  of  need 


Modern  rockets  operate  at  supercritical  pressures  with 
respect  to  the  propellants 

*  Need  to  understand  mixing  and  combustion  beyond 
liquid,  gas  states 

Shear  coaxial  injectors  are  a  common  choice  for  cryogenic 
liquid  rocket  engines 

Interactions  of  transverse  acoustics  with  injector’s  own 
modes  and  mixing  needs  to  be  understood  for  combustion 
instability 

*  Need  to  understand  differences  in  response  to  pressure 
and  velocity  nodes 

Understand  what  non-dimensional  numbers  capture  the 
mixing  of  typical  injectors 

Characterize  how  geometry  affects  mixing 
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Relevant  physics  of  shear  coaxial  jets 
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Transverse 
Acoustic  mode  from 
chamber/siren 

-  f=f(c,  geometry) 

Acoustic  modes  for 
outer  and  inner  jets 

-  f~c/2L  -  2  speeds  of 
sound  up  and 
downstream 

Wake  resulting  from 
inner  post  thickness 

-  St=fw/Uch 

Shear  layer 
instabilities 

-  Ste=f0/Uch 

Jet  preferred  modes 

-  St=fDij/Uij 
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Highlights  from  previous  work  on  jet 
instabilities  1/2 


Michalke,  1964 

—  Linear  stability  theory  for  inviscid  instability  of  a 
hyperbolic  tangent  velocity  profile 

Crow  and  Champagne,  1971 

—  Single  jet  preferred  mode,  Std=fd/U~0.3 
Ko  etal, 1976-1989 
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—  Some  of  earliest  detailed  description  of  near  field  mixing 
for  coaxial  jets 

Boldman  et  al,  1975 

—  Experimental  and  theoretical  analysis  for  mixing  of  two 
air  streams  with  different  velocities  -  points  out  different 
vortex  interactions,  St,~0.2  (Uave) 

Gutmarkand  Ho,  1983 

—  Collects  previous  results  on  jet  preferred  mode,  Std  has  a 
range  from  ~0.24-0.64 


Kwan  and  Ko,  J.  Sound  and  Vibration,  48  (2),  1976 


Wicker  and  Eaton,  AIAA  J,  (32)  No.3,1994 


•  Wicker  and  Eaton,  1994 


—  Forces  air  inner  and  outer  jets  independently  -  observes 
vortex  growth 

*  Dahm  et  al,  1992 


—  Seminal  pictures  of  different  instabilities  plus  effect  of 
absolute  velocity  and  R 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


Highlights  from  previous  work  2/2 


Balarac,  da  Silva,  Metais  et  al  (2003,  2007) 

—  DNS  analysis  of  coaxial  jets  -  same  density,  top-hat  profiles 

—  Consider  two  shear  layers,  study  effect  of  R 

^  ^  7777777771 

—  Consider  axisymmetric  and  azimuthal  excitation 
Buresti,  Talamelli,  Petagna  (1994, 1998) 

—  Air  jets,  same  density,  top-hat  profile,  Stdo=fdo/Uoj~0.3  to  1 , 
function  of  xlDt 

Segalini,  Talamelli,  et  al  (2006,  2011) 

—  Air  jets,  same  density,  top-hat  profile,  Stb(Mp)=f b( I i p)/Uaverage 
Birbaud,  Ducruix,  Durox,  Candel  (2006-2007) 

—  Single  air  jets,  low  Re,  laminar,  top-hat  profile,  subjected  to 
acoustic  modulation 

—  Systematic  study  of  effect  of  modulation  in  terms  of  Std,  St0 
Tshohas,  Canino,  Heister  (2004,  2009) 

—  2D  unsteady  CFD  for  LOX/H2  elements  but  non-reacting 

—  Unforced  behavior,  found  StMp=fd|OX/U|OX~0.1 0-0.25 
Richecoeur,  Scouflaire,  Ducruix,  Candel  (2006) 

—  Forced  transverse  acoustic  excitation  of  flames 


uj 


J77777777 


f  t  I  A  9 


Balarac  et  al,  Phys  of  Fluids  (19),  2007 


Segalini  et 
al,  Phys  of 
Fluids  (23), 
2011 


_ i _ i 

' _ 

1 

Birbaud  et 
al,  Phys  of 
Fluids  (19), 
2007 


Richecoeur  et  al,  JPP  (22)  No  4,  2006 
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Relevant  variables  for  cold-flow  studies 


Geometry 

Acoustics 

Recess 

Phase 

mm 

Coupling 

Single  Jet 

On/Off 

N/A 

v’  max 

N/A 

No 

LAR_thickLip 
(injector  1) 

1/2D1 

No 

SAR_thinLip 
(injector  2) 

1/2D1 

D1 

2-phase  P<Pc 

Off 

P’  max 

U’max 

0 

0.1-20 

SAR_thickLip 
(injector  3) 

0 

P>Pc  0.1-20 

T>Tc 

T<Tc 

Yes&  No 

LARJhinLip 
(injector  4) 

0 
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Experimental  setup  -  EC-4 
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Piezo-Siren 


Coaxial  Injector 


Inner 

Chamber 


Chamber  Pressure  vs.  Time 
(Subcritical,  J  =  9.6,  f  =3.01  kHz) 


Waveguide 


-Baseline 


Supercritical 
fluid  region 


Thermocouple  and  P 
transducer 


triple  point 
sublimation  line  J 


Temperature 
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Geometric  Rationale 


LAR_ThinLip 


r 


SAR_ThinLip 


LAR_ThickLip 


SAR_ThickLip 


LAR_thickLip 

SAR_thinLip 

SARthickLip 

LAR_thinLip 


D1 

(mm) 

D2 

(mm) 

D3 

(mm) 

D4 

(mm) 

t/D 

t=(D2- 

D1)/2 

Ao/Ai 

0.51 

1.59 

2.42 

3.18 

1.05 

12.9 

1.40 

1.65 

2.44 

3.94 

0.09 

1.6 

1.47 

3.96 

4.70 

6.35 

0.84 

2.9 

0.70 

0.89 

2.44 

3.94 

0.13 

10.6 

Two  recesses 


SAR,  LAR  ->  Small,  Large  Area  Ratio 
ThickLip,  ThinLip  ->  Post  lip  thickness 


LiJio  i  ■  oirticf  i  / 
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Chronological  progression  (only  coaxial 
results  are  summarized  in  what  follows) 
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*  Single  jets,  no  coaxial  flow 

—  Davis  et.  al.  (Ph.D.  thesis)  -  single  jets,  no  coaxial  flow 

*  Coaxial  jets 

—  Davis  et  al.  (Ph.D.  thesis)  -  LAR  thick 
—  Leyva  et.al.  -  LAR_thick 

—  Rodriguez  et.  al.  (Ph.D.  thesis)  -  LAR_thick,  SAR_thin 

—  Graham  et.  al.  -  SAR_thin,  two  recesses 

—  Teshome  et.  al.  (Ph.D.  thesis,  expected  April  2012)  - 
complete  all  four  geometries 

9  Also  complete  modal  analysis  of  earlier  geometries 

*  Future:  Combusting  coaxial  jets 

—  Wegener  et.  al.  (Ph.D.  thesis)  (in  process) 

—  Article  219  facility  funds 
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Dark-Core  Length  Measurement 
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•  First  raw  grayscale  images  were  converted  to  binary  images 
•  A  contour  was  drawn  around  the  “dark-column”  in  the  binary  image 
•  Axial  length  of  the  dark-column  measured  and  defined  as  the  Dark-Core  Length,  L 


Grayscale  (8-bit)  Binary  Image  Raw  Image  w / 
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Dark  Core  Length  -  SAR-thick  -  All 
data 
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Dark  Core  Length  without  Acoustics 


Outer  to  Inner  Jet  Momentum  Ratio  (J) 


■  Inj  1  Sub  Davis 

■  Inj  1  Super,  Davis 

♦  Inj  1  Sub,  Leyva 

♦  Inj  1  Near,  Leyva 

♦  Inj  1  Super,  Leyva 

▲  Inj  1  Sub,  Rodriguez 

▲  Inj  1  Near,  Rodriguez 

▲  Inj  1  Super,  Rodriguez 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


Dark  Core  Length  -  SAR-thick  -  All 
data 
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Dark  Core  Length  without  Acoustics 


■  Inj  1  Sub  Davis 

□  Inj  1  Tr<1  Near,  Davis 

■  Inj  1  Super,  Davis 

♦  Inj  1  Sub,  Leyva 

♦  Inj  1  Near,  Leyva 

o  Inj  1Tr<1  Near,  Leyva 

♦  Inj  1  Super,  Leyva 

A  Inj  1  Sub,  Rodriguez 

a  Inj  1  Near,  Rodriguez 

A  Inj  1  Tr<1  Near, 
Rodriguez 

A  Inj  1  Super,  Rodriguez 


Outer  to  Inner  Jet  Momentum  Ratio  (J) 
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Dark  Core  Length  -  SAR  thick 
supercritical  pressure 
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Dark  Core  Length  -  SAR  thick  - 
supercritical  pressure 
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Baseline  Dark-Core  Lengths 
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Acoustic  Field  Set-Up:  Pressure 
Antinode 


*  Pressure  antinode  (PAN)  -  condition  of  maximum  pressure  perturbation 
in  the  acoustic  field 

*  Piezo-sirens  forced  in-phase 

*  Superposition  of  quasi-ID  acoustic  waves  traveling  in  opposite 
directions  =>  PAN  at  the  jet  location  (geometric  center  of  test  section) 


t  (ms) 


t  (ms) 
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Image  interpretation  key 
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Baseline: 

Acoustics 

off 


0° 

Acoustics  on 
PAN,  VN 

“Pressure  coupled” 


180° 

PN,  VAN  " 
Largest  difference 

O 

expected  from  0 
“Velocity  coupled” 


► 


► 


► 


No 

dCC'J 

sties 


<T 


iatip 


pressure  =  fixed 
J  (p0u02/pjUj2)  =  fixed 

PN  -  pressure  node  -  Min 
PAN  -  pressure  antinode  -  Max 
VN  -  velocity  node 
VAN  -  velocity  antinode 
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Image  interpretation  key 
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pressure  =  fixed 
J  (p0u02/piU;2)  =  fixed 


Acoustics  on 
0° 

PAN,  VN 

“Pressure  coupled” 


Subsequent  images: 
acoustics  on; 

C 

Phase  increases  by  45 


Baseline: 

Acoustics 

off 


PN  -  pressure  node  -  Min 
PAN  -  pressure  antinode  -  Max 
VN  -  velocity  node 
VAN  -  velocity  antinode 
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180° 

PN,  VAN 
Largest  difference 

O 

expected  from  0 
“Velocity  coupled” 


360° 

should  look 
similar  to  0° 


-Critical  Pressure:  Two  Geometries 
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p 
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PN 
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J=p0U02/piUi2 
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ear-Critical  Pressure:  Two  geometries 

LARthickLip 


LAR thinLip 


Baseline 
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SARJhinLip 


Pr=1.05,  J=1.7,  p7p=0.32% 
Coupled  to  Outer  Jet  Acoustic  Mode 


*  Nearcritical  conditions 
*  No  jet  bending 
observed 
•  large  vortical 
structures 
generated  when 
coupled  to  injector 
mode 

*  Reduction  of  dark 
core  can  be  as  large 
as  90% 

•  More  clear  response 
of  jet  to  pressure 
antinode 


□ 


Subcritical  conditions 
Same  mode  - 
vortical  structures 
Not  as  dramatic 
reduction  as  with 
nearcritical  cases 


Pr=1 .05,  J=0.5,  R=2 
^^Pr=1.05,  J=0.5,  R=1 .3 


Pr=1.05,  J=2.2, 
p7p=0.25% 

Non  Coupled 
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Acoustic  Analysis  for  Injectors 


In  collaboration 
with  Jeff  Muss  and 
Rory  Davis,  Sierra 
Engineering 


Have  more  accurately 
computed  the  acoustic 
modes  for  the  inner  and 
outer  jets  for  constant  and 
linearly  varying 
temperatures  for 
subcritical  and 
supercritical  pressures 
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Amplitude  Amplitude  Amplitude 


New  analysis:  synchronized  p’  and 
images  taken  with  microscopic  lens 


1 

0.9  - 
0.8  - 
0.7 
0.6  - 
0.5  - 
0.4  - 
0.3  - 
0.2  - 
0.1  - 
0 


Chamber 
Center 


3000  6000  9000  12000 

Frequency  {Hz) 


3000  6000  9000  12000 

Frequency  (Hz) 


Inner  Jet  Pressure  Chamber  Wall  Center  Pressure 


Time  (ms) 


Outer  Jet  Pressure 


Time  (ms) 

Inner  Injector  Recessed  0.5  D. 

Pr  =  1.03,  J  =  5.1,  VR  =  5.9 
Forcing  Condition:  3000  Hz,  Maximum  Ap 


Uo=6.  25  m/s,  Ui=1.06  m/s,  To=165K  ,  Ti=118  K 
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Sample  Animation  -  PAN  (fF  =  3.14 
kHz) 


•  LARthin  Pr=  0.44, 
J  =  0.5 


Raw  Images 
Time  =  0.00  ms 


ROMs  1  and  2 
Time  =  0.00  ms 


I  I 
•  % 


POM  1 

Time  =  0.00  ms 


POM  2 

Time  =  0.00  ms 


t 

Superposition  of  POMs  1  and  2  Resulted  in  Downstream  Propagating 

Structures 

Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


CPSD  Magnitude 


Results  -  LAR-thin,  Pr  =  0.44,  Baseline  f  © 
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*  Antisymmetric  flow  structures  indicated  helical  type  flow  instabilities  for 

J=  2.1  J=  5.2  ,/=  11  J=  20 


Characteristic  peaks  broadened  and  shifted  to  higher 
frequencies  with  increasing  outer  jet  velocity 


■  oivtic*  I 
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CPSD  Magnitude 


Results  -  LAR-thin,  Pr  =  0.44,  PAN 


*  Gradual  shift  from  symmetric  to  antisymmetric  flow  structures  with 
increasing  J 

•  Response  at  forcing  frequency,  fF,  dominant  at  lower  J 

j=  2.1  J=  5.2  /=11  J=  20 
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CPSD  Magnitude 


Results  -  SAR,  Pr=  0.44,  Baseline 
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9  Helical  type  flow  instabilities  became  more  well-defined  with  increasing  J 


12 


4  D 


8  D 


J=  17 


J  =  2.0 


/(kHz) 


Unlike  LAR  flows,  characteristic  peaks  showed  minimal  variation  in 
_ frequency  with  outer  jet  velocity _ 
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CPSD  Magnitude 


Results  -  SAR,  Pr=  0.44,  PAN 


*  Symmetric  structures  persist  despite  increasing  J 


•  Response  at  fF  strong  at  highest  J 


J-  2.0 


J=  12 

(/f  =  2.92  kHz) 


J=  17 

(fF=  2.90  kHz) 


J  =  5.2 

(/f=  3.02  kHz) 


/(kHz) 


106 


/(kHz) 


/(kHz) 


Baseline 


PAN 
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Collaboration  with  ECP:  Thomas  Schmitt,  Juan  Rodriguez 
(AFRL  post-doc),  Sebastien  Candel,  Ivett  Leyva 


29 


20 


Longitudinal  slice  of  axial  velocity.  White:  minimum;  black:  maximum. 
Dark  line  indicates  iso-contour  of  zero  axial  velocity.  Top:  J  =  1.1;  middle:  J 

=  3.0;  bottom:  J  =  9.3. 


Top:  centerline  profile  of  density;  bottom: 
centerline  profile  of  velocity.  Dark  line:  J  = 
1.1;  dash  point  line:  J  =  3.0;  dashed  line:  J  = 
9.3. 


Velocity  Fields  for  Different  J  values 


600  = 


500  - 


j*  400 


^  300  - 


100" 
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Fundamental  frequencies  for 
baseline  conditions 


Case  “N6” 
(J=3.05) 
Density 
distribution 
(white:  60 
kg/m3; 
black:  410 
kg  /m3; 
logarithmic 
scale). 


OUTER 

JET 

J=1.05 

R=2.48 


Probe  A 


Case 

Ste  =hef/Uoj 
(probe  A) 

St'  =dif/Uij 
(probe  B) 

St'=l,f/UOJ 
(probe  C) 

N2 

0.25 

0.34 

0.15 

N6 

0.25 

0.26 

0.14 

Case  “N8” 

( J  =  9.3) 
Density 
distribution 
(white:  60  kg/ 
m3;  black: 
maximum ; 
logarithmic 
scale). 


OUTER 

JET 

J=3.05 

R=4.18 


Found  relevant  St  numbers  for  our  configurations 


I  ] 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


Combustion  Instability  Lab  - 
Background 


Combustion  instability  is  an  unsustainable  cimwth 
of  pressure  and_  heat  transfer  fluctuations  in  a 

rocket  engine 

—  Irreparable  damage  can  occur  in  <1s 

Combustion  Instability  caused  a  4-yr  delay  in  the 
development  of  the  F-1  engine  used  in  the  Apollo 
program 

—  More  than  $400M  for  the  propellants  alone  at 
2010  prices 

—  More  than  2000  full  scale  test 


Damaged  F-1  engine  injector 
faceplate  caused  by  combustion 
instabilities 


“Combustion  instabilities  have  been  observed  in  almost  every 
engine  development  effort,  including  even  the  most  recent 
development  programs”  -  current  JANNAF  Stability  Panel  Draft 
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*  The  single  largest  unknown  in  combustion  instability 
is  the  “combustion  response”  (how  combustion 
responds  to  acoustic  waves) 

*  Within  the  combustion  response,  flame  holding  in  the 
near  injector  field  is  a  key  mechanism 

*  GA  Tech  was  selected  to  lead  an  effort  on  a  closed- 
loop  study. 
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EC-4H  -  Combustion  Instability  Lab 


•  Concentrate  on  near  injector  field 


•  Measure  p’  and  q’  simultaneously  to 
evaluate  the  Rayleigh  criteria  for 
combustion  instability 

•  Start  with  shear  coaxial  jets  -  cold  flow 
heritage 


•  Lab  designed  for  2000  psi  -  about 
double  the  pressure  from  other  labs  in 
the  world 


•  Start  with  current  design  for  acoustic 
drivers  pressure  nodes  and  antinodes 
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•  1 500  g  LN2  tank  installed  Sept  28,  201 1 

•  Installed  Class  I  Div  2  outlets  (1 1 5V,  208V,  240V)  -  stripped  floor  for 
02  compatibility 

•  First  chamber  pieces  arrived  week  of  March  5,  201 2 

•  Fluid  systems  in  full  construction 


Construction/integration  to  be  finished  by  3Q  2012 
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*  Mixing  for  shear  coaxial  jets  has  other  major  variables  other  than 
momentum  flux  ratio  (J) 

—  Ratio  of  inner  jet  temperatures  is  another  important  variable 
—  Geometry  -  Area  ratio  and  lip  thickness  also  affect  mixing 

•  For  LARthick 


—  Found  bending  mode  with  largest  effect  at  velocity  antinodes 

—  The  reduction  on  the  dark  core  length  was  greatest  for  a  medium  J  range 

—  For  near  critical  pressures,  the  collaboration  with  ECP  determined  relevant  St  for  our 
injector  configuration  and  was  able  to  capture  qualitative  behavior  of  natural  and  excited 
jets 

•  ForSARthin 

—  Did  not  see  bending  mode  for  conditions  studied 

—  Saw  vortex  roll-up  and  puffing  occurring  over  entire  J  range  (0.09-21)  tested  - 
—  PAN  forcing  produced  symmetric  flow  structures  regardless  of  J 
—  Spectral  plots  showed  strong  response  to  PAN  forcing  at  low  and  high  J 

•  For  LARthin 


ui«  I  Cftwovin  ■  &tvu«c  I  (Hurt* 


—  Sees  both  bending  and  vortex  roll-up  modes  depending  on  the  acoustic  frequency 

—  PAN  forcing  at  low  J  produced  symmetric  flow  structures,  while  at  higher  J,  influence  of 
forcing  subsided 

I  —  Spectral  magnitude  plots  showed  decreasing  influence  of  PAN  forcing  with  increasi 
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AFRL  -  Edwards  was  tasked  to 
develop  an  experiment  and  to  lead  a 
joint  experimental  /  modeling  team 
to  study  coupled  flame  holding 
mechanisms. 

LEVERAGE  EXPERIENCE 
GAINED  WITH  COLD  FLOW 
EXPERIMENTS 
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•  First,  form  a  row  vector  consisting  of  all  pixel  intensity  values  of  each 
snapshot  image  (with  resolution  of  n  rows  by  m  columns)  in  order  of 
increasing  columns,  then  increasing  rows 


J  = 


1 

m 


n 

m 


i  =  l  /=1  i  =  2  /  =  2  i  =  n  i  =  n 

7  =  1  j  =  m  7  =  1  7  = m  7  =  1  7  =  ™ 


•  Then,  combine  all  such  row  vectors  for  N  sequences  of  image  frames 
resulting  in  a  matrix  A  consisting  of  N  rows  by  (M  =  nxm)  columns  of 

infpnsifv  valnec 
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Matrix  Decomposition 


•  Eigenvalue  decomposition  or  singular  value  decomposition  (SVD)  can  be  used 


•  SVD  Subroutine  readily  available  in  MATLAB® 


•  Prior  to  matrix  decomposition,  the  temporal  mean  of  A  was  subtracted  resulting  in  a  matrix  of  intensity 
fluctuations  A,  i.e.,  i 

a,=a.-t?Za» 

Ni  j=  1...M 


•  Application  of  SVD  on  A  gives  two  orthogonal  matrices  U  (NxN)  and  V  (MxM),  and  a  diagonal  matrix  S 
(NxM)  of  singular  values  in  increasing  order  of  magnitude 


A  =  USV  =  QVT 

•  Thus,  a  time-resolved  set  of  images  intensity  fluctuations  A(jc,t)  can  be  represented  as  a 
linear  combination  of  orthonormal  basis  functions  (j)k  such  that 

M 

A(x,t)  =  YJak(t)A(x) 

k- 1 


where  ak(t )  are  time  dependent  orthonormal  amplitude  coefficients  and  ^(x)  are  proper 
orthogonal  modes  of  A. 


•Equivlance:  columns  of  Q  ~  ak(t),  columns  of  V  ~  ^(x) 
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Acoustic  Analysis  for  Injectors 


In  collaboration 
with  Jeff  Muss  and 
Rory  Davis,  Sierra 
Engineering 


JUJSY5  11,0 

jah  12  2011 


Have  more  accurately 
computed  the  acoustic 
modes  for  the  inner  and 
outer  jets  for  constant  and 
linearly  varying 
temperatures  for 
subcritical  and 
supercritical  pressures 


LJJkCi  i  CHKOrvin  ■  SiVUCf  i 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


Acoustic  Analysis  for  Injectors 


In  collaboration  with  Jeff  Muss  and  Rory  Davis,  Sierra  Engineering 
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Have  more  accurately 
computed  the  acoustic 
modes  for  the  inner  and 
outer  jets  for  constant  and 
linearly  varying 
temperatures  for 
subcritical  and 
supercritical  pressures 
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Effect  of  R  on  L/D  for  a  given  J 
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Injector  IV:  Constant  J ,  Varying  R 


J  =2.0 

J  =2.0 

rrt 

II 

R  =  3.5 

L!D1  =  17.9 

1/D!  =  12.1 

J  =2.0 

J  =2.1 

• 

^r 

II 

R  =4.3 

L/D1  =  11.6 

L!Dl  =  10.1 
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Collaboration  with  ECP:  Grid  and  Mesh 


CENTRALE 


PARIS 
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Boundaries  with  modulated  normal  velocity 


3D  visualization  of  reservoir 


,=sin(2ir  f  t ) 


,  v2=sin(2irft  +  <p) 

Longitudinal  cut  of  domain  with  the  BC's 


Mesh  detail  near  the  injector 


CFD  Injector 


Solver:  AVBP  -  state  of  the  art  LES  code 

*  With  real  fluid  properties  to  tackle  supercritical  fluids 

Mesh:  2  100  000  nodes/10  000  000  tetrahedra 

*  Highly  refined  near  the  injector  (0.032  mm  on  a 
distance  of  10  inner  jet  diameters) 

CPU  hours  on  Europe  Supercomputer 
Center:  100,000 
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Effect  of  Temperature  on  L/D  for  a  given  J 


J  =  9.4  R=9.9 

109K/203K 

Ui=0.93/9.2m/s 

SARthinLip 


J  =  9.4  R=9.9 

128K/192K 

Ui=6.6  /36.2m/s 
LARthickLip 
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Flow  Structure  Comparison 
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Dark-Core  Length 
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Total  or  curved 
ength  (Lt)  = 

average  of  left 
and  right  sides 


P  =  1.50  MPa  ;  Pr=  0.44;  VR  =  7.50;  J  =  2.64;  acoustic  field  on 


•  Threshold  Images  based  on  Otsu’s  method  (  n.  otsu,  "A Threshold  selection 
Method  from  Gray-Level  Histograms, “  IEEE  Transactions  on  Systems,  Man,  and  Cybernetics,  vol. 
9,  no.  1,  pp.  62-66,  1979.) 

•  Accounts  for  variability  from  image  to  image  (including  dl  the 
parameter  by  which  the  jet  is  normalized) 
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Acoustic  Forcing  Results 
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Instantaneous  images  of  the  simulation  of  an  acoustic  case  with  the  injector  at  a  pressure  antinode  for  a  J  =  3.0 
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Dx=  0.70  mm 
D2  =  0.83  mm 
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D4  =  1.97  mm 
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Injector  exit  plane 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


Synchronized  p’  and  images 
Pr=1.03,  J=5.1,  VR=5.9 
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Synchronized  p’  and  images: 
Pr=1.03,  J=5.4,  VR=6.0 
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Sampling  Frequency=41kHz,  Driving  frequency=3.05kHz, 
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J  =  2.6  with  thick  inner  post  injector 
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Bending  vs.  vortical  structures  with  fine  atomization 
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^chamber  =  ^  MPa  (Pt  =  0.44) 


J  =  23  with  thick  inner  post  injector 


J  =  18  with  thin  inner  post  injector 


Two  different  mixing  mechanisms  for  similar  large  J  values 
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Crow  and  Champagne,  1971 

—  Single  jet  preferred  mode,  St=fd/U~0.3 
Ko  etal,  1976-1989 

—  Some  of  earliest  detailed  description  of  near  field  mixing  for  coaxial  jets 
Boldman  et  al,  1975 


Entrommeni  reqon 


Secorriory 

regioi 


—  Experimental  and  theoretical  analysis  for  mixing  of  two  streams  with 
different  velocities  -  points  out  different  vortex  interactions  Kwan  and  Ko>  J 


Sound  and  Vibration,  48  (2),  1! 


Gutmark  and  Ho,  1983 

—  Collects  previous  results  on  jet  preferred  mode,  St  has  a  range  from 
~0.24-0.64 
Dahm  et  al,  1992 


—  Points  out  importance  of  layer  thickness  and  velocity  defect  on  shear  layer 
seminal  pictures  of  different  instabilities 


Villermeaux  1998 

—  Inner  jet  core  length  expression,  LID^=6IJ0  5 
Richecoeur  et  al  (Candel’s  group),  2006 


Dahm  et  al,  JFM  (241)  1992 


—  Forced  transverse  acoustic  excitation  of  flames  | 


Richecoeur  et  al,  JPP  (22)  No  4,  2006 
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Subcritical  conditions 

•  Acoustics  undulate 

inner  jet  like  in 
Heidmann’s 
experiments 
consistent  with 
being  driven  by  an 
imposed  velocity 
field 

•  Reductions  of  core 
of  more  than  20% 

for  0.2  <  J  <10 
Near  and  supercritical 
conditions 

•  Less  undulation  on 
inner  jet  but  same 

mode 

•  Not  a  clear  effect  of 

phase  on  the 
reduction  of  the 
dark  core  length 
•  Max  reduction  of 
inner  jet  occurred  at 
1<J<5 


J  =  0.55  J  =  4.9 


VR  =  2.0  VR  =  5.9 
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Outer  Jet  Pressure 


Pr  =  1,03,  J  =  5,1,  VR  =  5.9 
Forcing  Condition:  3000  Hz,  Minimum  Ap 


Feb  18,  2010 
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'  attocube 


chamber 

+  Touterje,  (uPstream) 

+  T.  (upstream) 

inner  jet  ^  1  ' 

Outer  Injector  Post 
- Inner  Injector  Post 
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Effect  of  Recess:  SAR_ThinLip;  Pr=0.45,  J~0.09 


Flush  inner  post,  J=0.09,  p7p=0.45% 


Recessed  inner  post,  J=0.089,  p7p=0.60% 


Qualitatively  similar  at  very  low  J  values 
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